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1 Introduction 



The Rosenblatt process appears as limit of normalized sums of long-range dependent series (see [5], [2S]). In 
' the last years, this stochastic processes has been the object of several research papers (see [12], [17]j I1H]j 

among others). Its analysis is motivated by the fact that the Rosenblatt process is self-similar with stationary 
increments and moreover its exhibits long-range dependence (or long-memory). In this sense, it shares many 
properties with the more known fractional Brownian motion (fBm in the sequel) except the fact that the 
Rosenblatt process is not a Gaussian process. Recall that the fBm is the only Gaussian self-similar process 
with stationary increments. 

The practical aspects of Rosenblatt process are striking: it provides a new class of processes from which to 
model long memory, self-similarity, and Holder-regularity, allowing significant deviation from fBm and other 
Gaussian processes. The need of non-Gaussian self-similar processes in practice (for example in hydrology) is 
mentioned in the paper |26| based on the study of stochastic modeling for river- flow time series in |15| . 

The Hurst parameter H characterizes all the important properties of a Rosenblatt process, as seen above. 
Therefore, estimating H properly is an important problem in the analysis of this process. The Hurst parameter 
estimation from a A'^-length path of a fBm or more generally of self-similar or long-range dependent processes, 
has a long history. Several statistics have been introduced and studied to this end, such as parametric 
estimators (maximum likelihood estimator) as well as semi-parametric estimators (spectral, variogram or 
wavelets based estimators). Informations on these various approaches can be found in the books of Beran [3] 
or Doukhan et al. [10] . But in the particular case of the fBm, the main results concerning this estimation 
were certainly obtained by [TT] and [7]. Indeed the optimal method for estimating H is obtained from the 
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maximization of the Whittle approximation of the log-hkehhood introduced by Whittle in [30] (in the sequel, 
the Whittle estimator). This estimator shares the same asymptotic behavior than the maximum likelihood 
estimator (MLE), notably it is asymptotically efficient, but numerically the Whittle estimator is clearly many 
more interesting than the MLE (no need to inverse the covariance matrix). These properties also hold for long 
memory stationary Gaussian processes as it was established in and [7] under almost general conditions: 
the Whittle estimator is asymptotically normal with a ^fN convergence rate. 

In the case of long memory non-Gaussian time series, there exist very few results concerning the limit 
behavior of the Whittle estimator. In [TB], the case of fourth order moment linear processes have been 
considered and it was proved that the Whittle estimator is still asymptotically normal with a y/N convergence 
rate. But in [M], the cases of functionals of long memory Gaussian processes have been studied and the 
conclusion is different: in general, the Whittle estimator satisfies a non-central limit theorem with a non 
Gaussian limit distribution and a convergence rate smaller than \/N . This is notably the case when the 
functional is the Hermite polynomial H2{x) = a;^ — 1. 

Estimating the memory parameter of the Rosenblatt process appears to be a challenging problem. This is 
due to the fact that this process is not a Gaussian process, the explicit expression of its probability density is 
not know and standard techniques cannot be applied in this case. The development of new criteria, based on 
the Malliavin calculus and chaos expansion into multiple Wiener- Ito integrals (see the monograph ^\), for the 
convergence of sequences of random variables recently led to new results. In j23 and the authors studied 
the asymptotic behavior of the quadratic variations of the Rosenblatt process in order to obtain the asymptotic 
properties of an estimator for the self-similarity index. An approach based on wavelets has been also proposed 
in [3]. A common denominator of all these works is that the estimators constructed are consistent but in 
general their limit behavior is not Gaussian. This not very convenient for practical aspects. 

We want to put a new brick to the theory of the long-memory parameter estimation for non-Gaussian 
stochastic processes. We analyze the limit behavior of the Whittle estimator for the self-similarity index H of 
the increments of a Rosenblatt process (where 0.5 < H < 1). We will see that, as in the case of the estimators 
based on the quadratic variations (see [28]), the Whittle estimator has a non-Gaussian limit behavior and the 
convergence rate is N^~^ . This is due to the fact that, if one compares with the fBm case, the chaos expansion 
of the estimator involves a new term with a strong dependence structure that cannot be compensated by the 
smoothness of the estimator. This result is not totally surprising. Indeed, from [25], the Rosenblatt can be 
obtained as the limit of ^ X)fe!*i H2{Xk) when n — >■ oo, where {Xk) is a long memory Gaussian process and 
we know from [TJ] that the Whittle estimator of {H2{Xk))k also satisfies the same kind of non- central limit 
theorem. 

Unfortunately, our new result concerning the Whittle estimator keeps open the following question: is it 
possible to propose an asymptotically Gaussian estimator of the self-similarity parameter of a Rosenblatt pro- 
cess? However, Monte-Carlo experiments attests that the Whittle estimator numerically provides accurate 
estimations, clearly better than with other well known estimators and with almost the same quality than the 
one obtained for the fBm when H is close to 0.5. Hence, even if it asymptotically satisfies a non-central limit 
theorem, the Whittle estimator is really interesting for estimating the H parameter of a Rosenblatt process. 

We organized our paper as follows. Section 2 contains some preliminaries on multiple stochastic integrals, 
the Rosenblatt process and the Whittle estimator. In Section 3 we analyze the asymptotic behavior of the 
Whittle estimator for the self-similarity index. Finally, Section |4| contains a numerical study of the estimator 
and main proofs of this paper are established in Section [5] 

2 Preliminaries 

In this section we introduce the basic concepts used throughout the paper. We will present the the multiple 
stochastic integrals, the definition and the immediate properties of the Rosenblatt process and the basic facts 
concerning the Whittle estimator. 

To begin with we call back some elements about multiple stochastic integrals. 

2.1 Multiple stochastic integrals 

Let W = {Wt)t£R be a classical Wiener process on a standard probability space (Jl, J^, P). If / G L^(K") 
with n > I integer, we introduce the multiple Wiener-Ito integral of / with respect to W. The basic reference 
is the monograph [18]. Let / e Sn be an elementary function with n variables that can be written as 



/ = X^ii i '^ii,...,in^Ai^x...xAi^, where the coefficients satisfy Ci-^....^i^ = if two indexes ik and ii are equal 
and the sets S B{M.) are pairwise disjoint. For such a step function / we define 

where we put T4^(A) = /jj lA{s)dWs. It can be seen that the apphcation J„ constructed above from 5„ to 
L^(f2) is an isometry on 5„ in the sense 

E [In{f)Im{g)]=nl{f,g)L2(T^) ifm = n (1) 

and 

^ [/n(/)/rn(.g)] =0 if m 7^ n. 

Since the set 5„ is dense in L^(M") for every n > 1 the mapping /„ can be extended to an isometry from 

iv^(R") to L^(ri) and the above properties hold true for this extension. 

It also holds that /„(/) = /«(/), where / denotes the symmetrization of / defined by 



f{xi,...,Xn) = —^Yf{Xcr{l),...,X„(^n)), 



a running over all permutations of We will need the general formula for calculating products of 

Wiener chaos integrals of any orders m, n for any symmetric integrands / £ L^(M™) and g E L^(R"), which is 



mAn 

' m \ I n 



/rn(/)/45)= E^^L (2) 

1=0 ^ / V - / 

where the contraction f ®i g is defined by 

(/ ®l 9)isi, ■ ■ ■ , Sm-e,tl, ■ ■ ■ ,tn-i) 

/(si, . . . , Sm-e, Ml, • ■ ■ , ue)g(ti, . . . , tn-e, wi, ■ ■ ■ , ue,)dui . ..due (3) 

and by extension f ®o g = f ® g. Note that the contraction (/ (g)^ g) is an element of L^(M™+"~^^) but it is 
not necessarily symmetric. We will denote its symmetrization by {f®ig). 



2.2 The Rosenblatt process 

Recall that a fBm B = {Bt)t>o with Hurst parameter H £ (i, 1) and parameter C > is a centered Gaussian 
process with covariance function 

Gov (Bt, B,) = E (BtB,) ^ ^ {s^" + t''' -\t~.s\'"), s,te[0,oo), (4) 

with C = Var (Bi). It is the only normalized Gaussian _ff-self-similar process with stationary increments. The 
fBm admits the following moving average representation: for every t > and for every H > ^ 

Bt = Ccb{H) ^ " y)+'^dv^ dW{y). (5) 

where {W{y),y G M) is a standard Brownian motion with time interval R and cb{H) is a strictly positive 
constant that ensures that E{BtY = C^t^^ for every t > 0. 

The Rosenblatt process is related to the fractional Brownian motion. It shares many properties of the 
fBm. For instance, it has the same covariance function ^ as the fBm, it is _ff-self-similar and it has stationary 
increments. It has the same order of the Holder regularity of its sample paths as the fBm (that is, the 
Rosenblatt process, as well as the fBm, are Holder continuous of order 5 with < 5 < H). There are also 
some differences with respect to the fBm. One of them, is that it is defined only for the self-similarity index 
H > ^ and another difference, more important, is that it is not Gaussian. It can be expressed as a double 
multiple integral with respect to the Wiener process and therefore it is an element of the second Wiener chaos. 



More exactly, a Rosenblatt process Z — {Zt)t>o with self-similarity order H £ (i, 1) and parameter C > is 
defined by 

Zt = Ccz{H) [ I (I {u- yi)t~\u - y2)t'^dv) dW{yi)dW{y2) (6) 
where cz{H) is a strictly positive constant such that E{ZiY = C^, i-e. 

We will call a (standard) Rosenblatt random variable every random variable that has the same law as Zi with 
parameter C = f . In the sequel we will use the kernel Lh of the Rosenblatt process defined by 

L"{yi,y2)^ f {u~yi)f-\u-y2)f''du for (yi, ^2) G K^. (8) 







As the second order properties of a Rosenblatt process are the same than the ones of the corresponding 
fBm, the process Y — {Yt)t>o of the increments of a Rosenblatt process, defined by 

Yt = Zt+i - Zt for t > 

with Zq = by definition, is a long memory stationary process with covariogram 

rff,c(t):=Cov(yo,Ft) = ^(|t+lp^ + |t-lp^-2|tp^) for t e K, (9) 
and a spectral density fn.c defined for A e [— 7r,7r] by: 

f^^ciX) ^ E ^(^) = ^^^^(^fj^^"(-^) (1 - cos A) E |A + 2fc.r — , (10) 

fcez fcez 

since J^{1 - cosx)\x\-'^~^" dx = 27r(i7r(2ff) sin(7ri7)) (see Sinai, 1976, or Fox and Taqqu, 1986). 
2.3 The Whittle estimator 

Our purpose is to study the asymptotic properties of the Whittle estimator of parameters H and C computed 
from a sample {Zi,--- ,Zjv) of a Rosenblatt process Z — (Zt)t>o ■ Let us briefly introduce the Whittle 
estimator. The first step is to define the periodogram of the process Y of the increments of Z: 



1 1^-1 2 
^v(A) = ^|E^^^""' (11) 



fc=0 

2c 



Now, let 5 : R M be a 27r-periodic function such that g e L ([— tt, tt]) and define 

r ^ 

JN{g) = / g{X)InWdX, the integrated periodogram of y 

J —TV 

and J{g) ^ / g{X)fH,cWdX, 



with fn.c denotes the spectral density of Y defined in (fTO| . From ([9]) and (fTTj) . we also have 

^^W = ^ E ^Nik)e-^'' with r}v(fc) = - E YiYj+\k\. (12) 

|fc|<JV J = l 



which is a biased estimate of rH,c{k) (see the next section). Thus, the periodogram In{X) could be a natural 
estimator of the spectral density; unfortunately it is not a consistent estimator. However, once integrated with 
respect to some function, its behavior becomes quite smoother and can allow an estimation of the spectral 



density. The Whittle's contrast is a special case of the integrated periodograin JN{g)- Indeed, this contrast 
can be written (see [30] ) 



Un{H,C)^ r (log(/^,c(A)) + Jjv(l//if,c(A)))dA. 



The Whittle estimator is thus: 

{Hn,Cn) ■■= Aigmm(^H,c)eii/2,i)x{a,oo}UN{H,C). 

But using a classical renormalization (see Fox and Taqqu, 1986 for instance), the spectral density fH,c of Y 
can be decomposed as: 

t2 = ^ ^l!iyr(2i/)sin(7riJ) 

fH,cW = <y'9H{X) with<( ^''"^ (13) 

5h(A) = aff(l-cosA)^|A + 2A:7r|-i-2H 

feGZ 



dt 



for aU A e [-7r,0) U (0,7r] with a/f := exp -/ log HI - cosi) ^ |t + 2fc7r|"^~^^ j 
Then 

\og{gH{X))dX = (14) 

for H e (1/2, 1) and the minimization of Un{H, C) can be write again as a minimization in (H, cr^) and this 
implies that 



Hn = Argmin^g(i/2,i) {'^w(1/5h)} Argmin^g(i/2,i){y" ~~^'^^} 



Remark 1 However, for practical use, these definitions Hpf and a'j^ have to be modified. Indeed, the assump- 
tion that the process has zero mean is unrealistic. Moreover, the integrals defining the estimators has to be 
replaced by their approximations by a Riemann sum. Thus, define: 



= Argm^n^,^,/,^,^[- ^ A and <r^^-^. 



(vrfc/iV)/ '^"'^ N {^^gj,^{^k/Ny 



k=l 

With /at (A) - -—I {Yk ~Y N)e-'''^ for A e [-tt,^] and Yn^-Y.^^' 



N 

k=0 k=0 



3 Limit theorems for the Whittle estimator of parameters of a 
Rosenblatt process 

In this section, we prove limit theorems for the integrated periodogram JN{g) for 9 £ IL2[— 7r,7r]) satisfying 
certain conditions. Applied to g — l/gn, we will prove the almost-sure convergence of Hpf. Applied to 
g = -^[\/ gii) and then with a classical Taylor expansion also using the case g = -^r[j{^/ gn) this will provide 
a non-central limit theorem satisfied by Hn ■ 

Therefore, the main point is to obtain limit theorems for JN{g) ■ For this, we use a (now) standard approach 
based on the chaotic decomposition of the random variable Jn into a sum of multiple stochastic integrals. 
Since the Rosenblatt process at fixed time is a multiple integral of order 2, the product Y^lj+fc is a product of 
two double multiple integrals which can be expressed as a the sum of a multiple integral of order 4, a multiple 
integral of order 2 and a deterministic function, which is the expectation of YjYj^f^. This decomposition is 



transferred to the integrated periodogram J^r, which wih be written as sum of two multiple integrals (one of 
order 2, one of order 4) plus its expectation. What we do next, is to analyze these three terms that compose 
Jm- We will see that, as in the case of the variation statistic of the Rosenblatt process (see (25) or of wavelet 
statistic (see the dominant term is the one in the second Wiener chaos, which will give the asymptotic 
behavior of Jm- A detailed study of this term shows that it converges to a Rosenblatt random variable. 



3.1 Chaos decomposition of the integrated periodogram 

The purpose of this part is to provide the asymptotic behavior of a sequence Jn{9) using its chaos decompo- 
sition. We start first with the analysis of rN{k) with fixed fc G Z. We will study the convergence of rN{k) to 
the covariance function rn.cik)- We first observe that 

7V-|A;|-1 

and this converges to rH.c{k) when N — >■ cxi: rN^k) is an asymptotically unbiased estimator for rH,c{k). On 
the other hand, for every j > 1 we can write the increment of the Rosenblatt process as 

Y,^Z,+,-Z,=l2{ALf) 
where is the kernel of the Rosenblatt process (|8]) and we denoted by Aij^ the two-variables kernel 

Aif(yi,y2) = Ccz{H) {u~yi)f-\u-y2)f''du. 



In the sequel we will simply denote := L. Then, for every |fc| < N, by the product formula ^ 

N-\k\-l 

N 



rN{k)-rH.,c{k) ^ ^ J2 H^L,)h{AL^+\k\) 



W-|fe|-l 



^ ^ (/4(AL, ® AL,+|,|) + 4/2(AL, ®i AL,+|,.|)) 
+E?N{k)-rH,c{k). 



Therefore, for g e L^([— tt, tt]). 



, .TT N-\k\-l 

JM{g) -IEJn{9) - ^ / d^aW E E (^4(Ai, ® A,+|,|L) +4/2(AL, ®, AL,+|fe|)) (17) 

''^'^ \k\<N 1=0 

= T2,n{9) + T2m{9), 

1 r / X 

where T^^^^g) -.^ -— d\g{\) ^ e*'^'^ /4 (Ai,(yi, y2)AL,+|fc| (yi, 2/2)) (18) 

•'"^ \k\<N j=0 

2 r / r \ 

and T2,n{9):^^ d\g{X)Y, E e'""^ h[ dxl^L,{yi,x)l^L,+\k\{V2,x)). (19) 

\k\<N J=0 "'^ 

The notation T2^n{9), T4,n{9) suggests that the sequence T4jv(g) belongs to the fourth Wiener chaos while 
T2,N{g) belongs to the second Wiener chaos. 



3.2 Asymptotic behavior of the integrated periodogram 

In order to study the asymptotic behavior of Jn{9) let us first specify our assumptions on the function g. For 
H E (1/2, 1) and for g : [— vr, tt] i-> R, we introduce 



Assumption A{H): g E L^([— tt, tt]) and for any 6 > 0, there exists C such that 

\g{X)\ < C lAp^-i-'' for all A G [-tt, 0) U (0, tt]. 
A very useful point for us is the following (already proved in 12 , p. 531): 

Lemma 1 Let H e (1/2,1) and let gn be given by ilO\) . Then, l/gn, ^(l/Sff) '"^'^ '^ij^i^l 9h) satisfy 
Assumption A{H). 

The limit in distribution of the sequence JN{g)~E J Nig) as N oo depends on the asymptotic behavior of 
the two terms T2,N{g) and Ti^^ig) above. First, let us analyze the term in the second Wiener chaos, denoted by 
22,7v(<?)j of the decomposition Jjv(.9) when the function g satisfies Assumption A{H). Its asymptotic behavior 
is described by the following result. 

Proposition 1 For every H e (i, 1) and all function g satisfying Assumption A{H), 



with R a standard Rosenblatt random variable (IE R^ — 

Under the normalization of T2^N{g) the summand T^^^ig) goes to zero: 
Proposition 2 For every H E (i, 1) and all function g satisfying Assumption A{H), 

A non-central limit theorem satisfied by JN{g) is thus the consequence of both the previous Propositions 
□ andU 

Proposition 3 For every H E (^, 1) and all function g satisfying Assumption A{H), 



N'-" {Jn{9) - E JN{g)) 8 ^ ( .9(A)/(^+i)/2,c(A) ^a) R, 

with R a standard Rosenblatt random variable (IE R^ — 

The above Proposition [3] holds for any function g that satisfies Assumption A(H). The next result shows 
that, when g — (sh") (^^^^ 9h defined in (1131) ) under the normalization of T'2,7v(3) the deterministic term 
in the chaos expansion of JN{g) converges to zero. 



Proposition 4 With gn given by we have 



\dH^gH N~^oo 

It is also possible to show the following almost sure limit theorem for the sequence Jnig) when g satisfies 
Assumption A(H). 

Proposition 5 For every H E (^, 1) and all function g satisfying Assumption A{H), we have 



JNig) ^ a' / giX)gHiX)dX. 



Now, we can state our main result. 




Figure 1: FFT estimation (Silverman's method) of the density of the limit of Hn for H = 0.65, N = 5000 
from 1000 independent replications of Rosenblatt process paths. 

Theorem 1 Let be defined by ITS)) . Then 

Hn ^ H and N^-" (Hn - H) ^ -/{H) R 

TV-s-oo W-s-oo 

where R is a standard Rosenblatt random variable (with E R^ — 1) and j{H) is defined by: 

It is also possible to provide the strong convergence of Cn (defined in [TO]) to C: 
Corollary 1 Under the assumptions of Theorem]^ 

Cn ^ C and N^-" (Cn - C) A C p(H) R, 
where R is a standard Rosenblatt random variable (with IE R^ = 1) and p{H) is defined in (j42p . 

4 Monte- Carlo experiments 

We generated 200 paths of Rosenblatt processes for several values of H (= 0.55, 0.65, • • • , 0.95) and A'^ 
(= 1000, 5000 and 20000). These paths are obtained from the algorithm already used in [3j and deduced 
from the asymptotic behavior of ~ 1) when n ^ oo where (Yi) is a sequence of centered and 

normalized LRD processes with memory parameter d — H/2 (typically FARIMA(0, i?/2, 0) processes). Note 
that the generator of Rosenblatt process paths as well as the computations of the estimators used in this 
section are available on the websitelhttp: //samm.univ-parisl . fr/- Jean-Mar c-Bardet| with a free access on 
(in Matlab language). 

We computed the Whittle estimator Hn of H. An example of the estimation of the probability density function 
of Hn provided by the Silverman's nonparametric procedure is given in Figure [TJ This estimated probability 
density function appears could first appear as Gaussian density function but it is slightly asymmetric as the 
Rosenblatt density function should be when H — 0.65 (see for instance Figure 1 in |29j). 



The Whittle estimator Hn of H applied to paths of Rosenblatt processes is also compared to both other 
estimators: 

• The extended local Whittle estimator Hadg defined in [IJ from the seminal paper [53]. The trimming 
parameter is chosen as to = ]\j0.65 j^^j^jg jg an adaptive estimator) following the numerical recommen- 
dations of [T]. 



N = 1000 


H = 0.55 


H = 0.65 


H = 0.75 


H = 0.85 


H = 0.95 


mean Hjsr 


0.570 


0.653 


0.736 


0.815 


0.917 


std Hjv 


0.030 


0.041 


0.047 


0.053 


0.050 


mean Hadg 


0.570 


0.634 


0.708 


0.795 


0.906 


std Hadg 


0.072 


0.084 


0.094 


0.105 


0.102 


mean Hwa 


0.499 


0.542 


0.619 


0.685 


0.766 


std Hwa 


0.104 


0.116 


0.115 


0.129 


0.119 



N = 5000 


H = 0.55 


H = 0.65 


H = 0.75 


H = 0.85 


H = 0.95 


mean Hjv 
std Hjv 


0.582 
0.014 


0.655 
0.019 


0.743 
0.029 


0.837 
0.033 


0.929 
0.035 


mean Hadg 
std Hadg 


0.575 
0.041 


0.627 
0.052 


0.723 
0.062 


0.824 
0.067 


0.919 
0.072 


mean Hwa 
std Hwa 


0.550 
0.055 


0.610 
0.062 


0.698 
0.072 


0.800 
0.079 


0.891 
0.075 




N = 20000 


H = 0.55 


H = 0.65 


H = 0.75 


H = 0.85 


H = 0.95 


mean Hjsr 
std i^jv 


0.571 
0.008 


0.656 
0.015 


0.746 
0.020 


0.847 
0.025 


0.937 
0.025 


mean Hadg 
std Hadg 


0.563 
0.033 


0.637 
0.040 


0.734 
0.052 


0.838 
0.040 


0.931 
0.045 


mean Hy/a 
std Hwa 


0.569 
0.040 


0.630 
0.039 


0.728 
0.052 


0.838 
0.053 


0.931 
0.042 



Table 1: : Comparison between the Whittle estimator Hn and other famous semiparametric estimators of H (extended 
local Whittle estimator Hadg and wavelet estimator Hwa) applied to 100 paths of Rosenblatt processes with several 
H and values 

• The wavelet estimator Hwa as it was already defined in [3]. 

Note that both these estimators are semi-parametric estimators (while is a parametric estimator) and 
thus can be applied to other processes than Rosenblatt processes. However the asymptotic study of Hadg for 
Rosenblatt was not already done. The results are detailed in Tabled] 

Conclusions of the simulations: The results obtained with the Whittle estimator Hn are very convinc- 
ing and fit the limit theorem we established. Indeed the estimator Hn seems to be asymptotically unbiased 
and its standard deviation is depending on H (the larger H the larger the standard deviation). 
For giving a comparison, it could be interesting to compare these results with those obtained for fractional 
Brownian motions (fBms) instead of Rosenblatt processes (the computation of the Whittle estimator is exactly 
the same for both those processes). In additional simulations, we obtained the following results: the standard 
deviation sn of Hn applied to fBms almost not depends on H. Hence for almost all H e]0.5, 1[ and = 1000, 
SN ^ 0.02, for N = 5000, sat ~ 0.01 and for N = 20000, sn ^ 0.005. As a consequence, from Table [H 

• when H — 0.55, the standard deviation of Hn for Rosenblatt process is close to the one obtained for 
fBm. In terms of theoretical convergence rate this is not surprising since the convergence rate of Hn is 
^0.45 jjj^ (,g^gg q£ Rosenblatt process while it is N'~'-^ in case of fBm. 

• when H = 0.95, the standard deviation of Hn for Rosenblatt process is dramatically larger than the one 
obtained for fBm. Once again, this is not surprising since the convergence rate of Hn is N'^-'^^ in case of 
Rosenblatt process while it is still iV°-^ in case of fBm. 

Also from Table [1] we can compare the convergence rates of Hn and both the semiparametric estimators 
Hadg and Hwa- The standard deviation of Hn is almost the half of those of Hadg and Hwa- The results 
are really convincing and show the good accuracy of the Whittle estimator for estimating H . However, to be 
fair, we have to underline that Hn is a parametric estimator typically appropriated to Rosenblatt (or fBm) 
increment processes while Hadg and Hwa are semiparametric estimators which can be applied to general 
classes of long-memory processes. 

Remark 2 As we improved the generator of Rosenblatt processes, the results obtained with Hwa o-i^e o, little 
better than those obtained in especially when H is close to 1. 



5 Proofs 

The following technical lemma will be needed several times in the sequel. 

Lemma 2 Let g satisfy Assumption A{H) with H G (1/2, 1). Then, for any S > and € Z, 

5(A) e^'^^dX = 0((1 + \k\)-^"+^). 
Proof: This have been stated and proved in [12], Lemma 5. 

Proof of Proposition [1} First, using the definition of the contraction (E)i (see (jSj), 
(ALj (g)i ALj+|fe|)(yi,y2) = dxALj{yi,x)ALj+m{y2,x) 



C Cz{H) dx du{u-yi)l {u-x)l I dv{v-y2)l {v - x) 

Jr Jj J j+M 

C^cl{H)fi{l-H,-) duj^^^^ dv{u-yi)l \v~y2)l \u~vf-^ 



where we changed the order of integration and we used the identity: for (a, h) E (—1, 0)^ and a + b < —1, for 
(m, v) e with V, 

/uAv 
[u ~ yfiv - yfdy = /3(-l - a - 6, a + 1) - vr+'+\ (21) 
-OO 

Then 



T2Mg)-C \ ^ ^ > dXgiX) J2 E ^'"'M du dviu-y^)J~\v-y2)^'\u-v 

By Lemma[3l the sequence {N-^^^T2.N{g))N has the same limit in L^(51) as the sequence {N^~^T2 j^{g))is! 
where 

N-\k\-l 

ikX 



fi»(.)=c'|e|="^/'.A,(A)x: j: 

PU 2 > \k\<N j=0 

, rJ + 1 rJ + \k\+l rj+l „ „ 

x/2(/ du dv\u^v\"-' du'{u'-yi)l \u'^y2)l ). 



But from the definition of a Rosenblatt process and the computation of /j*^^ du j^^'!^^^^^ dv\u — v\^ ^, 



?2,^(5) = Cy^^^-L ^ (|%(A)e»'^^dA)(|fc + iri + |A-iri-2|fcri)i? 



where {Rt)t is a standard Rosenblatt process with parameter H. But from Lemma [21 with 5 > i) that can 
be chosen arbitrary smaU, we have /^^ 5(A)e*'=^dA = 0((1 + |fc|)"^^+'^') and since (|fc + + \k- - 

2|fc|^+^) ^ i(i/ + 1)H k^~^ when fc is large enough, we deduce that 

E ( r ffWe^'^A) + +\k- - 2ifci^+i) < OO. 

Moreover, since a Rosenblatt process is a i?-self-similar process, (iV^^i?Ar_|fc|)|fc|<Ar ^ (-Ri_ifci)|fc|<Af- This 
is also a continuous process and therefore 



^.H^, V ,^2 2(2i/-l) 1 



(.g) 



^""" uli+Hy h ^ (/%(AK'='rfA)(ifc+ir+^ + ifc-ir^-2|fci-+i)7?^ 



But (|fc + + |fc - - 2\k\"+'^) = 2r(jf+i)/2,i(fc) (see ©) and therefore from the definition of a 

spectral density, 



5(A)/(H+i)Ac(A)dA i? 



Lemma 3 Define the sequence of functions {Gn)n by 

N-\k\-l 

GM{yi.V2):=N-" / g{X)d\ ^ e^^^ 



Then, 



k\<N j=0 



GAr(j/i,2/2) 0. 

TV— i-oo 



J + 1 



(u'-yi)l 'K- 2/2)1 'du' 



Proof: 1. On the one hand, define the sequence of functions {Gn.i)n by 

Af-|fc|-l 



\ki\<N |fe2|<iV' 
N-\k^\-lN-\k2\-l -^ + 1 rn + \ki\ + l 



Gtv4(2/i,2/2) := f g{X)dX ^ 6"*^'= ^ 

/ du dv\u-v\ {u-yi)^ {v~-y2)l -{u-yz), 

Jj J] + \k\ L 

Then we have 

\\GNAyi,y2)\\l2(^.)=N-^" I dyidy, ^ ^ T .g(Ai)e-^^^^^dAi T g(A2)e-^'^'==dA2 

dvi\ui - vi\"-'^{ui - - y2)+^^ - (mi - 2/2)+"^]) 

Ff-i — -ir — -1 — 

dv2\u2-V2\ {u2-yi)l [{V2~y2)l - {u2 - Vl) I Jj- 

Then, using (HI]), there exist Di> Q and 1)2 > such as. 



Jl=0 J2=0 •'-'I 

X I / dM2 

'j2 -'i2 + |fe2 



ji+\kl\ 
j2 + |fe2|+l 



,-iA2fe2^\ 



< 



_D2 



E 



JV-|A;2|-1 jj+i + + i + 1 ^i2 + |fe2l + l 

X ^ Y, / / '^'"2 / |W1 - Wi|^"^|m2 - 1'2|^"^|W1 - W2I 

jl=0 ^2=0 ^1 + "'i2 "'j2 + |fc2| 

X [\U, - 1^21^-1 - |«1 - V2\"'' - \U2 - Vi\"-' + \VI - V2f~'] 
N 



H-1 



1 



|felMfc2|<W 



[(l + |fci|)(l + |A:2| 



]2H-S 



/ duidvidu2dv2 (|wi - wi - |u2 - t'2 - IfeH |ui - "2 + ji - ^2!) 



\ui -U2+ jl - j2\ - \ui -V2+ jl - j2 - Ifell - \U2 - I'l + j2 - jl - \kl\\ + \vi - V2 + jl - 32 + |fci| - | fc; 

^iJb E E Y^^i^^M) 



H-l 



\H-1 



(22) 



\ki\<N \k2\<N \j\<N 



with 



H-l 



I{ki,k2,j) ■■= r., , I, IW-, , I, |M2g-A / duidvidu2dv2{\ui -Vi - \ki\\ \U2 -V2 - \k2\\ \ui - U2+j\) 
[(1 + |fci|)(l + \k2\)]^" J[o^i]i 

\U1 - U2+j\"-^ - \UI -V2+J- \k2\\"-^ - \U2 -VI- J- ^f'^ + \vi - V2 + J + \ki\ - Ifcall^^^ 

(23) 

and where both the last inequahties are obtamed from changes of variables, Lemma [2] (with 5 > Q which can be 
chosen arbitrary small) and symmetry properties. For \x\ > 2 and (u, v) £ [0, 

and for < 2 and {u, v) S [0, 1]^, u ^ v, \x + u — v\^~^ < \u — v\^~^. Moreover from a usual Taylor expansion 
of the function x i-^ \1 + x\^^^ , there exists C > such as 



\y + j\"-'-\3\"-\l + {H-l)^) <Cj-^ when2|y|<j. 
Therefore, there exists C > such as 

\UI -U2+ ]\"-' - |U1 -V2+J- \k2\\"-^ - \U2 -VI- J- \ki\\"-^ + \vi - V2 + ] + \ki\ - I^H^-^ 

<C{l+\k,\ + \k2\)\j\"-' 

<c[{i + ifcii)2«-i+^- + (1 + ifc2i)^^-^+'] \jr"-' 

when |j| > 2max(|fci|, |fc2|). Then for |j| > 2max(|fci|, |fc2|) and for \ki \ and |/c2| large enough. 
Therefore, 

]^ E E E nkuk2,j) ^-^^ 0. (24) 

\ki\<N \k2\<N |j|>2max(|fci|,|fc2|) 

Now, for 2 < \j\ < 2max(|A:i|, |/c2|) and for |fci| and \k2\ large enough, 

X {\3\"-' + \UI -V2+J- \k2\\"-' + \U2 -V,-J- \k,\\"-' + |«i -V2+J + \k,\- \k2\\"-') . 

But for l3 e R*, since 1/2 < F < 1, there exists C > such as for < M, 

M ^ 

Thus, for 2 < |j| < 2max(|fci|, |fc2|) and for \ki \ and |/c2| large enough. 



y 



[(l + |fc,|)(l + |fc2pi^+^-^ (niax(|fci|,|fc2|))^^-^ 



Therefore, 

]^ E E E (25) 

\ki\<N \ k2\<N 2<|j|<2max(|fci|,|fe2|) 

We can easily add to this asymptotic behavior the cases j = 0, j = 1 or j = — 1 thanks to the convergence of 
the integral defined in [0, 1]^. Finally, we obtain: 

||Giv,i(yi,y2)lli.(K2) 0. (26) 



2. On the other hand, define the sequence of functions (Gjv,2)jv by 



JV-|fc|-l 



GnAvi^V^) ■■= r 9i\)dX J2 E 



|fc|<Af 3=0 

du I dv I du'\u-v\"^' {u-yi)l "(^-^2)+ yi)+ '("'-2/2)+ 

Jj+|fe| L 

Then following the same kind of computations and expansions than for \\GN.i{yii J/2)||^2(r2')j we have 

Ar-|fei|-17V-|/c2|-l „^ 



fel|<Af |fc2|<Af Jl=0 J2=0 
+ 1 rjl + \ki\+l rji+l I'j2 + l I'j2 + \k2\+l rh+l 

X / dui i dvi I du'i / du2 / dv2 / c?M2 jui — t;i|^~"'^|u2 — W2|^~"'^ 

•^jl -'il + |fel| Jjl J32 Jj2 + \k2\ J 32 



\2H-2 I / / |2H-2 

'■I ^ "2I 



^ ]^ E^^^ + |fc^|)(i + |fc,|)]2H^. y^d^id.id<dz.2rf^;2rf4(iui - + i^iii -V2 + \kr"-' 

I „,/ |2_ff-2 



\ki\,\k2\<N 



\3\<N 

< ^ V V V - (i + l?!)'''"' 

\ki\<N \k2\<N \j\<N^^ ' ^'■'^ ' "^'^J 

As a consequence, since all the previous sums are finite, we deduce that 

\\GN.2iyi,y2)fmR2) o. (27) 

Then ^ and jST]) imply Lemma [3] holds. ■ 

Proof of Proposition [2l In a first step we use the isometry of multiple integrals (H]) and the fact that 
the norm of the symmetrized function is less than the norm of the function itself. In a second step we 
change the order of integration and we use ([2T|). Then, for the last bound below, we use the same lines as in 
the proof of Lemma [31 We have 



IE (N^-"f4,N{g)f < N-^" I dyidy^dy^dy^ f dXid\2g{\i)g{X2) E E e-^'^^^^e-^'^^^^ 



JV-l/cil-l Ar-|fc2|-l +1 ^.ii+lfeil + l 

■r-^ ■r-^ /-Jl+J- Wl-l-l'-ll+J- _H__j^ H__^ H__^ H__-^ 

/ / dv{u-yi)^ {u-y2)+ {v~y?.)+ (w - 2/4)+ 

ji=0 32=0 "'ji "'ji + |fci| 

+ 1 (■32 + \k2\ + l „_ _H__ 

du' dv'{u'~y,)l \u'~y2)l [v - ys)^ {v'~yi)l 

1 32 J 32 + \k2\ 

P\l - H, |) N-'^ E E ( r ^''''''9i^i)dXi) ( r e-^-^^g(A2)dA2) 

|fei|<A'|fe2l<JV "'-^ "'-'^ 

Ar_|/ci|-iAr-|fe2|-i ji+i /•i2+|fc2|+i 

^ E E / / / / dMdvdu'd«'|M-uf^->-«f^-2 



Jl=0 J2=0 



< ^N'-'" E E „,,2H-. E k«aO>«a(j + |fci|-|fe|)|. 

|fci|<Ar|fe2|<w ((1 + Fil)(l + F2I)) IjKat 

As in the proof of Lemma [3l we decompose the previous right left hand term in two parts; firstly, when 
|j| > 2max(|fci|, \k2\), we have thAj) < c(l+ and thAj + M - Ifel) < c{l + \3\Y"-^. Then, with 



c> 0, 

N'-'" E E „ „ ,,,2H-S E \rHAj)rHAj + \ki\~\k2\)\ 

\ki\<N \k2\<N (,(1 + iklDi,!- + 1^21)) 2max(|fei|,|fe2|)<b"|<JV 

<cN'-'" y y jw^N*"-^ 

\ki\<N \k2\<N 

<cN^"-^, (28) 

since 2H — S > because H > 1/2 and S > can be chosen arbitrary small. 
Secondly, when |j| < max(|A;i|, \k2\), using \rH,i{j + |fci| — |fc2|)| < 1, we obtain 

N'-'" E E „ „ ,,,2H-s E \rHAj)rHAj + \ki\^\k2\)\ 

\ki\<N \k2\<N i(l + \kl\){l + 1^21)) |j|<2max(|fci|,|fc2|) 

<ciVi-2« V V 1 __ nrax(|fci|, 1^21)2^-1 

\ki\<N \k2\<N 

<cN^-^"+^. (29) 
As a consequence, from p8| and p9)) . 

JE (7Vl-^f4,Ar(5))' 0, 
\ / N-^oo 

and the conclusion of Proposition [2] follows. ■ 

Proof of Proposition [4) The proof follows the lines of end of the proof of Theorem 2 in [12], p. 528- 
529. Notice that 

3 ^ r 



-»(^(^))^E/_/^^(^)^^».c(^)» 

Denote 
Then 

From (3.3) in [12], we know that 



yrH^cik)ek=0. (30) 



feez 

Therefore, 



|fc|<iV 

= TVi"^ ^ rff,c(fc)efe-iV-^ ^ krH,c{k)ek. (31) 

fc|<Af |fc|<Af 

From Lemma[5]and since rn.cik) = 0(|fcp^~^), |fc| — > oo, we have for all (5 > 

rff,c(fc)efc = 0(|fcr2+^). (32) 

Then, from a usual comparison with a Riemann integral, X]|fc|<iv krH,c{k) = 0{N^). Moreover, from pop . 
S|fe|<Af "i^H.cik) Cfe = — X]|fc|>Af '''H,c{k) Sk and also from a usual comparison with a Riemann integral and ([5^ . 
we know that '}2\k\>N ^H,c{k) et — 0{N^^^^). As a consequence, from ([3T]) . 



As S can be chosen arbitrary small, we deduce that Proposition |4] holds. 

Proof of Proposition [Sj Note that the Fourier coefficients of In are given by 

e''''^lN[x)dx = fjv(fc)l(|fc|<jv)- 



The proofs of Proposition [T] and [5] imply that 

TMik) "^^^ rn.cik). (33) 
The above convergence holds almost surely by an argument in 9 (Theorem 7.1, p. 493). H 
Proof of Theorem [Jl For establishing the strong convergence, define for h £ (1/2, 1), 

VN{h)^^JN{l/9h). (34) 

t2 



From Proposition [SI VN{h) ^^—4 V{h) = a'^ J__^ gniX)/ gh{X)dX. As H — Argmin,jg(j/2,i)^(^) ^^id H, 
Argmin^jg^ ;2 ]^)V^r(/i), from usual arguments (see for instance |12 ), then Hn H. 



For proving the non-central theorem, apply Proposition |3] to g = (gTS") ^^'^ from Proposition |4l we obtain: 



with R a standard Rosenblatt random variable. Now use the following result established in Lemma 2 in |12j . 
Suppose that {bN)N is a sequence of real numbers such that — > + oo. Assume that 

N—^oo 

bN^VN{H) A Y, (35) 
dh Af^oo 

where K is a random variable. Then 



Using bN = " ,Y ^ d, \J H("-H) ' i ( S% f(H+i)/2..cW Ri, Proposition [3] applied to ^(^) and 
Proposition 01 then (1551) holds. This implies Theorem [1] ■ 

Proof of Corollary [1] Using VAr(/i) defined in ([M)) . it is sufficient to write 

Cn = {^l{HN)VN{HN)f'\ 

with 

^^^^-'""hrmk^) f-'^^(iAi), (36) 

and the strong consistencies established in Proposition [5] and Theorem [1] allow to show the almost sure con- 
vergence of Cm ■ 

For establishing the non-central limit theorem, the Taylor's formula implies that 

with probability 1, and with — H\ < \Hn — H\. From previous Theorem [TJ it follows 

7Vi-^(Ujv(i/) - CT^) = N'-"{Vn{Hn) ~ cT^) + Op{N"-^). (37) 



Moreover, from Proposition [3l 

N'-"{VN{H)-]E{VNiH))) A -C'(3{H)R. (38) 

with /3{H) = ^ H%"Hr ( I-.9H'Wf{H+i}/2,iWdx)- Since (Fjvli?)) = +o(iV-V2) (g^^ instance 
[H]), from dnZl) and dSHl), we deduce: 

7Vi-^(VA,(i?Ar)-a2) A i/3(iJ)i?. (39) 

From Theorem [T] and using the Deha- Method, 

N^-"{n{H^)~ti{H)) A (M'(i/))S(i/)i?. (40) 

Thus, using (l39l) and (l4Qj) . we obtain: 

iVi-^(Mi?iv)VA,(iIjv)-M^)^') C2((^'(i/))2^(i/)-i^ + l/3(iJ)Mi?))i?. (41) 

As a consequence, using again the Delta-Method, we have: 

N'-"{dN~C) -^^(^(^'(H)r^iH) + ^|3{H)^^^H))R, (42) 

achieving the proof of Corollary [1] 

■ 
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